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1. INTRODUCTION 

Perturbative QCD calculations give a good de- 
scription of data on jet observables as far as a 
scale of few GeV's is involved in the process. For 
instance, in e~^e~ annihilation the dependence of 
the mean jet multiplicity on the resolution pa- 
rameter yc defined in the Durham algorithm is 
well reproduced by QCD calculations in absolute 
normalizations down to yc ^ 10~^, i.e. to a scale 
of a few GeV's, when leading and next-to-leading 
terms in In yc are resummed and a matching with 
the full two- loop result is performed [Q. In the 
softer region with scales of the order of 1 GeV 
or less, deviations from purely perturbative pre- 
dictions are visible and this is usually explained 
in terms of purely non perturbative hadroniza- 
tion effects. However, following the idea of a soft 
hadronization mechanism a simple picture 
of hadronization, known as local parton hadron 
duality has been indeed formulated, where 
the same perturbative description of QCD parton 
cascade is stretched down to small scales of a few 
hundred MeV for the transverse momentum cut- 
off Qq, and then parton predictions are directly 
compared to hadron distribution up to an overall 
normalization factor. The interesting feature of 
this picture is that it is surprisingly successful in 
describing the phenomenology of hadron produc- 
tion, as far as sufficiently inclusive observables, 
like particle multiplicity or inclusive energy spec- 
tra, are concerned |^,^. 



Even though both parton multiplicities and jet 
multiplicities are described with the same type 
of evolution equation, two different sets of pa- 
rameters, and in particular, two different normal- 
ization factors, were needed so far to describe 
phenomenologically the experimental results on 
the two observables. Here we show that both 
observables can be described in an unified way, 
if the complete QCD master equation which de- 
scribes the evolution of a parton cascade is ex- 
actly solved. We refer to for further details 
on this calculation. 

2. THEORETICAL FRAMEWORK 

The perturbative QCD picture gives a proba- 
bilistic description of the production process of a 
multiparton final state in terms of a parton cas- 
cade. Formally, the physical information about 
this process is entirely contained in the evolution 
equation for the generating functional of the mul- 
tiparton final state 1^; in its more general for- 
mulation, this equation correctly takes into ac- 
count angular ordering, energy conservation and 
the running coupling at the one-loop order; the 
choice of the thansvere momentum as the vari- 
able which enters in the coupling incorporate at 
least part of higher order effects. The infrared 
cutoff k± > Qo also provides one with a natural 
regularization of the infrared singularity in the 
running coupling. By appropriate differentiation 
of the equation for the generating functional, a 



coupled system of two evolution equations for the 
mean parton multiplicities Mq and Mg in quark 
and gluon jets can be derived |,|l|,§: 



dz 



drj ~" 27r '-2 

{Ugiv + lnz)+Ug{7j + ln(l - z)) - Ugiv)} 

{Ngif] + lnz)+Afg{r^ + ln(l - z)) - J^giv)}] 
dMM f'-'^ «.(fc±) ^ . . 

{Mgiv + lnz)+Mg{v + ln(l - z)) - AA,(ry)}. 

The $AB are the AP splitting functions for the 
process A ^ B, Nc and n/ are the numbers 
of colours and flavours respectively, 77 — In ^ , 

where k = Q sin ^ j with Q the hard scale of 
the process {^/s in e+e^ annihilation) and 8 
the maximum angle between the outgoing par- 
tons. as{k±) is the one- loop running coupling, 
given by as(k±) = 2Tr/{b\n{k±/A)) with b ~ 
{llNc — 2nf)/3 and a smooth matching of the 
running through the heavy quark thresholds 
The argument of as is chosen to be the "Durham" 
transverse momentum 



k± = min(z, 1 — z)k, 



(2) 



< i, one 



The boundaries of the integral over the par- 
ton momentum fractions z are determined by 
the lower cutoff k± > Qc and given by Zc = 
(Qc%/2)/g = V2y^ = e-''. Since 
finds ?/c < f and then 77 > In 2. 

The coupled system of evolution equations 
must of course be supplemented by a suitable set 
of initial conditions; they are fixed at threshold 
to be: 



Afg{r]) = J\fg{r]) = 1 for < < In 2. 



(3) 



This master equation contains the correct lead- 
ing and next-to-leading terms in an y^al expan- 
sion, but it has been obtained by using approx- 
imations valid in the limit of soft emission (for 
instance, the definition of k± given in (||) comes 
from the Durham definition of k 1 in the limit 



of small 2.). In other words, non-logarithmic 
terms have been neglected in the derivation of 
the master equation, and these terms give a non- 
negligible contribution in the region of small rj 
(i.e., poor resolution (large yc) or low cms en- 
ergy (small Q)). To improve the description in 
this region, the contribution of 0{as) in (1) has 
been replaced by the explicit result for e^e^ — > 3 
partons in ©(a^ )[§§]: 

(4) 

where A/'q^'' is obtained by taking the first itera- 
tion of (1) with the initial conditions (^) and the 
full 0{as) contribution is given by: 



A^'''"*(yc) = 2 dz, dz2Q{d2z - yc) 

Jl/2 Jl-zi 



CFa.s{kA_) zf + zj 

2tt (1-Zl)(l-Z2) 



(5) 



^23 = min ( — , — ) (1 - zi) > j/c (6) 

where zi (22) denote the quark (antiquark) and 
z^ — 2 — zi — Z2 the gluon momentum fractions 
[Cf - 4/3). 

3. PHENOMENOLOGICAL ANALYSIS 

3.1. Multiplicities in e+e annihilation 

The average jet multiplicity at Q =91 GeV as a 
function of the resolution parameter yc — Q^c/Q^^ 
defined via the Durham algorithm, is shown in 
Figure 1|^ (lower data points). Also shown in 
the lower part of the Figure is the theoretical pre- 
diction in absolute normalization of the jet mul- 
tiplicity at the same cms energies of 91 GeV, as 
obtained from (^) . The predictions depend on the 
single parameter A only, so we can use the LEP-1 
data on jet multiplicity to fix the best value of the 
QCD scale A. The best value of A turns out to 
be 500 ± 50 MeV, in qualitative agreement with 
previous results within the MLLA framework ||l|. 
The main improvement of our new approach in 
comparison to the old results is that we can now 
describe experimental data down to very small 
values of the resolution parameter yc] therefore. 
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Figure 1. Data on the average jet multiplicity at Q = 91 GeV and the average hadron multiplicity 
(assuming N = ^Afch) at different cms energies with Qc = 0.507 GeV as a function of j/c- The solid lines 
show the predictions for the hadron multiplicity (upper curve) and for the jet multiplicity at different 
LEP cms energies (91, 133, 161 and 172 GeV) (lower curves), obtained by using eq. (^) with parameters 
Kail = 1, A = 0.5 GeV and A = 0.015. The right most data point for hadrons refers to pions only. 



after having properly taken into account the ef- 
fects of energy conservation, the validity of the 
perturbative picture can be extended well below 
the scale of 2-3 GeV which was usually claimed to 
be the soft scale where hadronization effects be- 
yond perturbative QCD take place. In the Figure 
the predictions on the jet multiplicity at LEP- 
2 cms energies are also shown. These curves 
could be compared with new data recently col- 
lected at CERN. After having fixed the first pa- 
rameter A, we can now ask ourselves whether it 
is possible to describe also the cms energy de- 
pendence of hadron multiplicity within the same 
framework. In this case, we can still choose two 
free parameters, the hadronic scale related to the 
k± cutoff Qo and the normalization parameter 
Kail, which relates parton and hadron multiplic- 
ities as Mall = KaiiMcorr ■ I* tums out that a 
common description is indeed possible, with pa- 
rameters given by Kaii = 1 and A = In Qq /A = 



0.015 ±0.005 (i.e., Qo 507 MeV). Figure 1 (up- 
per data points) shows the data on hadron mul- 
tiplicities from y/s = 1.6 (only pions) up to 172 
GeV p taken as JVaii — \Nch as a function of 
the resolution parameter yc = Qq/Q^, together 
with the perturbative predictions of (Q) with the 
aforementioned parameters. It is clear from the 
Figure that the perturbative approach can de- 
scribe quantitatively the hadron multiplicity even 
at very small cms energies, thus extending again 
the range of validity of the perturbative picture. 
We refer to for a more detailed comparison 
of our new numerical results with previous ana- 
lytical approximate solutions. Let us stress that 
a common description of both hadron and jet 
multiplicities is possible with the common nor- 
malization Kail = 1, which is a natural value, 
since it provides the correct boundary conditions 
at threshold for both hadrons and jets. 

It is interesting to study the dependence of our 
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Figure 2. Same data as in Fig. 1; The solid lines show predictions for the hadron multiplicity and for the 
jet multiplicity at 91 GeV with the same parameters as in Fig. 1; the dashed line show predictions for 
the parton multiplicity corresponding to the same A and Kaii but different values of the cutoff parameter 
Qo= 2, 1, 0.75, 0.66, 0.55 GeV (from bottom to top). 



results on the infrared cutoff Qq] one could indeed 
stop the parton evolution at a larger scale of the 
order of 1-2 GcV and then fill the gap between 
partonic predictions and experimental data using 
purely nonperturbative hadronization models, as 
implemented in the most used Monte Carlo mod- 
els. The perturbative predictions for the parton 
multiplicity with different values of the infrared 
cutoff Qq arc shown in Fig. 2. Notice that in our 
picture most of the particles (~ 3/4) are produced 
in the very last stage of the evolution for 0.5iQoil 
GeV, where Monte Carlo models already stopped 
the perturbative phase. A different value of Qq 
also changes the energy dependence of the pre- 
dicted parton multiplicity. It is indeed remark- 
able that the right energy dependence of particle 
multiplicity can be reproduced with a small value 
of the infrared cutoff Qo and a normalization fac- 
tor Kail = 1- A residual correlation among the 
two parameters exists, so a variation of one of the 



two parameters within 30% can be accomodated 
by adjusting the other one. 

It is important to stress that the evolution vari- 
able ?7 = In ^ in (1) is given by the ratio of the 
two dimensional scales which enter in the prob- 
lem, i.e., the external scale Q and the resolution 
scale Qc- Therefore, the same evolution equa- 
tion controls both the cms energy dependence 
of the cluster multiplicity, with clusters defined 
at a given resolution (where Q varies and Qc is 
fixed), and the resolution dependence of the jet 
multiplicity at a fixed hard scale Q (where Q is 
fixed and Qc varies). There is indeed only one 
place where the two scales enter separately and 
not throTigh their ratio, namely the expression 
of the running coupling ofg, where only Q (and 
A) but not Qc appears. Therefore, if one would 
switch off the running of the coupling, the equa- 
tion would depend only on the dimensionless vari- 
able r] and there would be no difference at all be- 
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Figure 3. Same as in Fig 1, but with A = 0.35 GeV and a different definition of k± (see text for details). 



tween jet and hadron multiplicities. It is then 
clear that the difference between the two multi- 
plicities entirely comes from the running of the 
couphng. Notice that this difference goes up to 
a factor 10 in the region of large yc (i.e., small 
cms energy for hadron multiplicity), and this fac- 
tor is within this picture completely related to 
the increased coupling at small transverse mo- 
mentum. This result further confirms previous 
results on energy spectra||,|l^. The effect of the 
running of the coupling is visible not only in the 
difference between the two curves, in particular 
at large values of Uc, but it is also directly visi- 
ble in the strong rise of the jet multiplicity in the 
high resolution region (low yc) (see Fig. 2). In- 
deed, the jet multiplicity increases roughly by a 
factor 3 if one lowers the resolution cutoff from 0.9 
down to 0.5 GeV, as expected by the quick grow- 
ing of the running coupling close to the Landau 
singularity (which is however screened by the Qc 
cutoff). Such results are qualitatively consistent 
with expectations from the Double Log Approxi- 
mation!^. 



3.2. Ratio of hadron multiplicities in quark 
and gluon jets 

The numerical solution of the coupled system 
of evolution equations (1) provides one with new 
results on the mean multiplicity Mg in gluon jets, 
and then also for the ratio of hadron multiplicities 
in quark and gluon jets, r = Ng/Nq. This observ- 
able has received a lot of attention both theoreti- 
cally and experimentally. Theoretically, different 
predictions were derived for the mean parton mul- 
tiplicity in a full hemisphere of the gluon jet pro- 
duced by a primary gg colour singlet state. The 
asymptotic value of r = 9/4 has been successively 
reduced to smaller values, by including next-to- 
leading order and even part of higher order cor- 
rection (see 1^,0 for a discussion of this point). 
Experimentally, a completely inclusive configura- 
tion similar to the theoretical predictions has re- 
cently been realized in the decay Y — > gg^ by the 
CLEG GoUaboration and in e^e" — > 3 jets 
with nearly parallel q and q recoiling against the 
gluon by the OPAL Gollaboration[|lO|. 

By numerically solving the evolution equation 



(1) with initial condition (without any finite 
order correction term yet), the prediction for the 
ratio r is further reduced with respect to the pre- 
vious approximate calculations and can describe 
the OPAL data for jets of 40 GeV very well. One 
fails, however, to reproduce the CLEO data for 
jets of about 5 GeV. This failure could be due to 
an important contribution from the fixed order 
correction term, not yet included so far. 

3.3. Another definition of transverse mo- 
mentum 

All results we have shown so far have been ob- 
tained by using the Durham-inspired definition 
of transverse momentum given in (j^). With this 
method, a value of the parameter A consistent 
with the previous finding from the energy mo- 
ments |8| has been found; however, the scale A 
(or, equivalently, Qo) has turned out to be larger 
than the value of ~ 250 MeV used in the descrip- 
tion of energy spectra|H,^,D. In the latter case, 
the alternative definition of fcj_ = z{l — z)k has 
been used. To check whether this difference plays 
an important role, an alternative calculation with 
this different definition of the transverse momen- 
tum has been performed. As shown in Fig. 3, 
a good description of experimental data is indeed 
possible also in this case, but with a lowered scale 
A ~ 0.35 GeV and the same A parameter. The 
remaining difference to the earlier A comes from 
having taken the new scale k = ^/2E instead of 
K = E in the new calculation. Therefore, the new 
results presented here are completely consistent 
with previous studies of the energy spectra. 

4. CONCLUSIONS 

The complete QCD evolution equations for 
quark and gluon jets has been numerically solved. 
By including also the full 0{as) correction for 
e^e~ annihilation, a common description of the 
resolution dependence of jet multiplicity at LEP- 
1 and of the cms energy dependence of hadron 
multiplicity in the whole cms energy range has 
been achieved with two free parameters only and 
an overall normalization factor consistent with 
1 . The ratio of hadron multiplicity in gluon and 
quark jets at LEP-1 is described as well. For such 



inclusive observables the perturbative description 
turned out to be valid also in extreme kinematical 
domains, thus extending the region of validity of 
the perturbative approach supplemented by the 
Local Parton Hadron Duality picture. Further 
tests of this picture and, in particular, of its va- 
lidity for soft particles and for less inclusive ob- 
servables are certainly needed; some of them are 
discussed in 
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